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' Abstract. Let be a bounded convex domain in Euclidean n space, 

X £ dfl, and r > 0. Let u = {u^ , . . . , u^) be a weak solution to 

V-{\Vuf~Wu) = in nriB{x,4r) with IVul*'"^ u„ = on dnnB{x,4r). 

We show that sub solution type arguments for certain uniformly ellip- 
tic systems can be used to deduce that | V{t| is bounded in n B{x, r) 



with constants depending only on n, p, N. and finB(x r)| ' argument 
replaces an argument based on level sets in recent important work of 



< 

■ [CM], [CMl], [GS], [M], [Ml], involving similar problems 
> I 

^ 2000 Mathematics Subject Classification. Primary 35J25, 35J70. 

Keywords and phrases: p harmonic systems, Neumann problem, p Lapla- 
cian, Gradient bounds, convex domains. 

> 

l> , 

O , 1. Introduction 

O 

ly-^ I Let X = {xi, . . . ,Xn) denote points in Euclidean n space, M", and let 

■ (•, •) denote the standard inner product on M". let \x\ = denote 
the norm of x and set B{x,p) = {y : \y — x\ < p} when p > 0. Given an 
open set O C M", let C^{0) denote infinitely differentiable functions with 
compact support in O. If 1 < g < cx), let W^''^{0) denote the Sobolev space 

^ I of functions g with distributional derivatives gxi,^ ^ i ^ n, and norm, 



\s\\w^,i{o) = \\9\\li(o) + ll|Vg|||L9(o) < 



oo. 



Here Vg denotes the gradient of g and || • || l<?(o) is the usual Lebesque q 
norm relative to O. li E,F C let H{E, F) denote the Hausdorff distance 
between the sets E,F, and let l-El denote the Lebesgue n measure of E 
whenever E is measurable. Throughout this paper we assume that 0, C M" 
is a bounded convex domain and for given x ^1,1 < p < oo,r > 0, that 
u = : fi— )']R^ is a weak solution to the p Laplace systems 

equation in f2 n i3(x,4r) with vanishing Neumann data on (?f2 n i?(x,4r). 
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That is, u'' G W^'P{nnB{x,4:r)),l < k < N, and if cj) = {(p^ , . . . 
(1.0) / \Vu\P-^ {Vu, V(t>)dx = 

JnnB{x,4r) 

whenever ^* e W^'P{Q fl B{x, 4r)) and 1 < i < iV, vanishes outside a set 
whose closure is a compact subset of B{x,4r). In the above display, Vu is 
the nN tuple, (Vtt^, . . . , V-u^) while the inner product is relative to R"^. 
In this note we show that 

Theorem 1.1. Let u,x,r,Cl,p be as above. There exists C depending only 
on n,p, N, r"'/\il fl B(x, r)|, such that for every y E Cin B{x, r), 

\Vu\P{y) < Cr-"" ! \Vu\Pdx 

JnnB{x,4r) 

We note that in [CM], [CMl], the authors studied weak solutions to quasi- 
linear elliptic equations and systems of the form V- (a(|Vtt|)V?i) = f{x) in a 
convex domain Q under both Dirichlet and Neumann boundary conditions: 
n = 0, = respectively on dQ. It is not obvious to us that their nearly 
endpoint global type results imply Theorem 1.1. Our proof though uses 
some of the same arguments as in these papers, including fundamental use 
of an inequality in [G] for convex domains. Also as in these papers, we first 
prove Theorem 1.1 for a weak solution to a related PDE in a smooth convex 
domain and then use a limiting argument to get Theorem 1.1. However the 
above authors use the level sets of \Vu\ and perform some rather involved 
calculations on these sets in order to obtain their results. In contrast we use 
well known sub solution type arguments for \Vu\ (i.e, Moser iteration) to 
get Theorem 1.1. 

Acknowledgment: We would like to thank Prof. N Garofalo for bringing 
the reference [CMl] to our attention and showing interest in our work. 



2. Proof of Theorem 1.1 

To begin the proof of Theorem 1.1 let c be a positive constant, not neces- 
sarily the same at each occurrence, which unless otherwise stated may only 
depend on n,p,N. We note that the p Laplace equation is invariant under 
dilations, translations, and scaling. Thus it suffices to prove Theorem 1.1 
when 

(2.0) r = l,x = 0, and / \Vu\p dx = 1, 

JnnB{o,4:) 

as we can then transfer back to the general case using the transformation 
y = X + rx, after multiplying u by an appropriate constant. We continue 
under assumption (2.0). We may also assume that O.ridB{0, t) ^% whenever 
t < A since otherwise it follows from (1.0) with ^ = a suitable extension of 
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u to i?(0,4) that u = constant. Let p, 3 < p < 4, be such that 

(2.1) / \Vu\Pda<c 

JnndB{o,p) 

where a denotes HausdorfF n — 1 measure = surface area. Existence of 

p,c follows from (2.0), writing the integral in polar coordinates, and the 
usual weak type estimates. Moreover, (2.1) holds for p belonging to a set of 
positive measure in [3,4]. Let r]{x) = 7(|x|) where 7 = — 1 in [0,/9 — (5], linear 
in [p — 5, p] and 'j{t) = for all t > p. Using this 77 as a test function in (1.0) 
and letting 5 — > 0, we deduce from the Lebesgue Differentiation theorem, 
that we may assume 

(2.2) / \Vu\P-^Uuda = 

JhndB{Q,p) 

where ^{x) denotes the outer unit normal to O fl B{0, p) at x eCih dB{0, p) 
and Uy = (n^, . . . ,u^) with u'l{x) = (^u^{x),v{x)), 1 < k < N. Next given 
e > 0, < e < 1/2, let = il(e) be a convex domain with 

(2.3) ncn, H{dn,dn)<€, and dnec°°. 

Let p be as in (2.2) and let / = (/\ . . . , /^) be defined by f{x) = 0, 1 < 
k < N, when x € H dB{0,p) while / = |Vtt|^~^tti/, otherwise on 

Q n dB{0, p). Note from (2.2) that / satisfies the compatibility condition : 

(2.4) f fda = 0,1 <k<N. 

Using (2.4) along with well known trace theorems and variational methods, 
we deduce the existence of a unique u = u{-, e) : $7 fl B{0, p)^R^ with 
g l^i'P(i7 n B{0,p)), l<k<N, satisfying ^^^(o.p) ^^^^^ = and 

(2.5) / {e + \Vu\y^^-\Vu,Vv)dx = [ {f,v)dx 

whenever v = {v\...,v^) with v'' € W^^f{nnB{0,p)). In (2.5) we have 
used V OIL ^Ci dB{0,p) to denote the trace of v. More specifically, u is the 
minimum of the functional 

Fe = - [ (e + \Dw\^y/^dx - [ <f,w> 

P JnnB{o,p) JdB{Q,p)nn 

among all w € W^'^ such that /j^n^^g p ^ ~ ^- Existence follows from lower 
semicontinuity of F^, (2.4), and compactness of the trace operator. Choosing 
V = u in (2.5) and using (2.4), Poincare's inequality, (2.1), we deduce that 

(2.6) / (e + I Vu|2)f 1 Vtx|2 dx < c. 
Jb{q,p) 

Next we claim that 

(2.7) v!" e C°°{nnB{0,2)) for l < A; < at and l < m < n. 
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Constants in (2.7) however may depend on e,p,n,N and the smoothness 
of dQ in (2.3). We sketch the proof of (2.7) in the appendix to this paper 
for the reader's convenience since we have not been able to find a suitable 
reference. Our proof uses a reflection type argument, after straightening d^l 
in an appropriate way, to first get C^'" regularity in the closure of J7nB(0, 2). 
After that we use Schauder type arguments to bootstrap, as in [ADN] , 
[ADNl]. 

From (2.5) and (2.7) we see that n is a strong solution to 
(2.8) V • ((e + |V«|2)f/2-iVu) = 

in Q f] B(0, p), where V- denotes the divergence operator. Let A G with 
|A| = 1. Let ux = (Vu^, . . . , Vuj^) denote the directional derivative of u 
in the direction A. Differentiating (2.8) with respect to A we get for fixed 
1,1 <l <N, that is a solution inQn B{0, p) to 

(2-9) EE^(^^r4(<))=o 

m=l i,j=i * 

where 

(2.10) b^{x) = (6 + I Vup)^/2-2 ((p - 2)n- < + S^^ie + \Vuf)) (x) 

when X E fin B{0,p). Here Sf^ = 1 if i = j,l = m, and is zero otherwise. 
We note that if C = i^k) is an n x matrix, |^|^ = (^,^) where the inner 
product is relative to M"^ and h = {e + |Vup)V2^ ^^^^^ aX x enn 5(0, p) 

N n 

(2.11) min(p-l, 1)1^1^/1^-2 < ^ <max(p-l,l)|ep/i^-^ 

m,l=l i,j=l 

For 1 < z < j < n let 

AT 

(2.12) Cijix) = (e + \Vu\\x)r'[{p - 2)(^ ^ .) + Sij{e + |Vnp) ] 
where 6ij is the Kronecker 6 and let L be the differential operator, 




We observe that if /x = (/xi, . . . , /x^) G M", then 

n 

(2.13) min(p — 1, l)|/xp < Cjj(x) /Ltj/Ltj < max(p — 1, 1) |/Lt|2 . 

We use (2.12), (2.13) to give a proof of a well known sub solution inequal- 
ity for L(/t^). To this end let be the point in with a one in the k th 
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position and zeroes elsewhere. Ifg>pata;Gr2n B{0,p), we calculate 

n 

(2.14) 

n 

Moreover, 

n N 

(2.15) L{hP)=P E E(^^'"'««-.)-^+^(^'') = ^i + ^(^'') 

i,i,k=l m=l 

where A is the Laplacian. Thus 
(2.16) 

n AT 



i,j,k=l m,l=l 



n N 
i,j,k=l m,l=l 

Using (2.9) with k playing the role of j and xj the role of A we have 

n N n N 

(2.17) ^4 = -pEE<.(/^'"'^^.-,)-^ = -Hhn+pY. E^'^'K-.)'- 

1=1 1=1 i,j=l 1=1 

Combining (2.15)-(2.17) it follows that L{hP) > p min(l,p - l)/iP-2|V/ip 
and thereupon from (2.13), (2.14) that for some c = c(p) > 1, 

(2.18) L(/i«) > c-^g2/i«-2|V/ip = (4c)-^|V(/i«/^)p. 

Let (p G C^(5(0,3/2)) and put w = (j)'^. From (2.18), (2.3), (2.7), and 
integration by parts we find for q > p a constant c = c(p, n) > 1 with 
(2.19) 

T =c-^ [ 02|V/i«/2|2dx < / (t>'^L{hi)dx 

if2nB(o,4) 7nnB(o,4) 

= E / <^'cy (/i9),,.z/i da - 2 V / #x,Ci,(/i«)^^dx = Ti + Ts. 

Prom (2.5), smoothness of dO- and it, we deduce that tijy = on 5r2nS(0, 2). 
Using this deduction and (2.12), we see that only the 5ij term contributes 
to the sum defining Ti. This remark and [G]( pp. 133-137 and Equation 
3.1.1.8) yield 
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„ n N 



(2.20) 

N 

where M(., ) is the second fundamental quadratic form on dO. and Vfu is the 
tangential component of V^x on d^. Since J7 is convex, M(., ) is nonpositive. 
Thus Ti < . Also using Cauchy's inequality with (5's we get 

(2.21) \T2\ <sj <^2|V(/i«/2)pdx + \V4>\''hidx 

Choosing S so small that the first term on the righthand side is <T/2 and 
then using (2.21), (2.20), in (2.19) we conclude after some arithmetic that 



(2.22) J (t)'^\Vh"l'^\^dx <cj |V</.p 



h'idx 



where c = c{p,n) > 1. From (2.22) and Sobolev's inequality for fl 5(0,2) 
appUed to 6' = h^/^cl) we deduce for n / 2 that 

(2.23) 

<(c/|J^nB(0,2)P)||V^||i.(^nB(0,2)) 

<{c/\i^nB{0,2)f) [ \V(l)fh'idx 
JnnB{o,2) 

If 71 = 2 replace 2n/{n — 2) by 4 in (2.23). We can now use Moser iteration 
(see[GT, ch 8]) in a well known way and get for every x G $7 fl B{0, 1) that 

(2.24) (e + |Vnp)f/2(x) <C [ (e + \Vu\y/'^dx 

where C has the same dependence as in Theorem 1.1. Finally, we note 
(see [D], [L], T] for = 1 and [Tl] for > i) that Vu(-,e) is Holder 
continuous on compact subsets of D B{0,p) with Holder constants inde- 
pendent of e. Also from (2.6) wc find that {u{-,€)} is uniformly bounded in 
VF^'^(J7 n B{0,p)). Using these facts it follows easily that subsequences of 
{u{-, e)}, {Vu{-, e)} converge uniformly on compact subsets of fl B{0, p) as 
e^O to u' € W^''P{^ n B(0,p)). From uniform integrability type arguments 
we see that (2.5) holds for u' , Q D B(0, p) with e = when v is infinitely dif- 
ferentiable on R". Since a given v G PF^'^(0 n 5(0, p)) can be approximated 
arbitrarily closely in the norm of this space by such functions, we conclude 
that (2.5) is valid with e = and u, Q. replaced by u' , Now u is the unique 
function (up to a constant) having these properties so tt' — u = constant. 
Thus (2.6), (2.24) hold with e = when u, Q are replaced by tt, fi. Hence 
Theorem 1.1 is true. □ 
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3. Appendix: Proof of 2.7 

After a rotation if necessary, around a neighborhood of any boundary 
point yo G dU D B{0,2), dO, can be locally represented as {y',yn) '■ Un = 
c^{y')} and Q = {(y',y„) : yn > Hv')} for some G As in [GT] 

(Section 14.6), it follows that there exists /i > such that x & Q and 
d{x,dQ) = d{x) < fj,, then there is a unique point y{x) € dfl, such that 
\y — x\ = d{x, dil) = d{x). The points x and y are related by 



(3.1) x = y- v{y)d 

where v is the outer unit normal to O 

As proved in Lemma 14.16 in [GT], the map g{y',d) = {y',4>{y')) — 

i'{y' ,(l){y'))d is locally invcrtiblc in a neighborhood U of {yQ,0) and maps 
?7 n {d > 0} into Hence by the inverse function theorem y and d in (3.1) 
are locally C°° functions of x. Using Vd{x) = —i^{y) one now calculates 
that the inverse of g is the mapping, 
(3.2) 

g ^(x) = (y' , d) where d = d{x) and yi = Xi — d{x)dxi{x), 1 < i < n — 1. 

By shrinking U if necessary, we may assume that g{U ri{d > 0}) is contained 
in Qn 5(0, 5/2). 

We note from (3.2) that at x G g(U), 

(3.3) {Vyi, Vd) = {si - {Vd^^)d - d^^Vd, Vd) = -d{Vd^„Vd) = 

where in the last inequality we have used the fact that Vd is constant along 
the line from x to y in (3.1) (so the directional derivative of dx^ in the 

direction of Vd{x) = 0). 

Define u{y',d) = u{x) where u is as in (2.5). We have 

n— 1 n— 1 

(3.4) {Vu, Vv) = Uy,Vyi + uaVd, K^Vi + ^d^^) 

i=l i=l 
Using the orthogonality of Vy^ and Vd, we get 

n-l 

(3.5) {Vu, Vv) = UdVd + ^ Uy^Vy.qij{y', d) 

In conclusion, by renaming the variable d = yn and using the change of 

variables formula, we deduce from (2.8), (3.5) that u satisfies 

(3.6) 

/ (e+ < A{y)Vu, Vu >f'^-^ < A{y)Vu, VC > V{y', yn)dy = 

where U is as above and C, is in Wq'^\U). V is the Jacobian of g. The matrix 
A = {aij) is defined by ajj = qij when 1 < i,j < n — 1, where qij is as 
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in (3.5), akn = o^nk = when k < n and Unn = 1- Also if z G M", then 
A{y)z = rj where 77 = (ryi, . . . , 77^) with 



r]i = ^ aij{y)zj and AVu = {AVv}, . . . , 

Now let ^ = (V. Using ^ in place of ( in (3.6), we obtain 
(3.7) 

(e+ < A{y)\7u, Vn >)p/'^-^ < A{y)Vu, > dy 



L 



Ur\{yn>0} 

N 

= (e+ < A{y)Vu, Vu >)p/2-i V < A{y)Vu^, VV > i^V'^dy 

Now we extend A, u and V to U r\ {yn < 0} by even reflection and denote 
them by B, v and W respectively. That is, B{y', yn) = A{y', -yn),v{y', yn) = 
Hy', -yn), and W{y',yn) = V{y' , -yn) for y.„ < 0. We note that B and W 
are Lipschitz extensions of A and V respectively and v € W^''p(U) 

Now for any C € Wl''P{U) 
(3.8) 

/ (e+ < B{y)Vv, Vv >)f/2-i < B{y)Vv, VC> dy = [ + I 

Ju Jun{y„>o} Jun{yn<o} 

For the first integral on the right hand side of (3.8), we use (3.7) and the 
fact that V, W restricted to {yn > 0} equal u, V respectively. For the second 
integral, we define ^(y',— 2/n) = Ciy'iVn) and change variables. Using the 
definition of v, i/j and W we obtain (3.7) with ^ replaced by tp. Altogether 
we see that 
(3.9) 

(e+ < B(y)Vv,Wv >)p/2-i < B{y)Vv,VC > dy 



L 



L 



u 

N 



(e+ < B{y)Vv, >)p/^-^ ^ < B{y)Vv\ VW > C}W~^dy 



^ 1=1 
Thus V G W^'P{U) is a weak solution in [/ to a system of the form 

(3.10) V-S(y,Vi))-So(y,Vi)) = 0. 

Moreover this system satisfies the structural assumption in [Tl] (see (1.1) 
and (1.7-(1.13) in this paper) except that in [Tl], the analogue of b is as- 
sumed to be in y. However estimates in [Tl] only use Lipschitz norms in 
y so are also valid in our case. To check the lower order term observe that 
the I th component of bo{y, Vv), 1 < I < N, is given by 

(3.11) (e+ < B{y)Vv, Vv >)p/2-i < B{y)Vv\\7W > W~\ 
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This term may be discontinuous in y when j/n = thanks to Wy^, but still 
satisfies the growth and structure assumptions in (1.13) of [Tl]. Therefore, 
we conclude from [Tl], that -0 G C^^'^iU) fl Wf^'^iU) (for second derivative 
estimates, sec section 4 in [Tl]). 

We note from the above results and (3.10), (3.11), that bo{y,'Vv{y)) is 
Holder continuous in U since Vy^ = at points in U where y„ = 0.. Moreover 
{bij)yf. are Lipschitz continuous when k < n. Therefore, by using a difference 
quotient argument we deduce that if A = j/jt for < n then for 1 < m < A^, 
(3.10) can be differentiated with respect to A as in the derivation of (2.9) 
in order to obtain that = (Vi)"^, X),l < m < N, is a weak solution to a 
uniformly elliptic linear system in N equations of the form, 

N n n 
m=l i,j=l i=l 

for fixed I, where fl, C^™', 1 < i,j < n,l < l,m < N, are Holder continuous. 
We can now apply Theorem 2.2 in Chapter 3 of [Gi] to conclude that vx 
has Holder continuous derivatives in U , i.e Vy^yj is Holder continuous in U 
when i + j < 2n. Now coupled with the fact that v is in W^''^{U) we can 
write (3.10) in non-divergence form and obtain for each I = 1, N that 



TV 



(3.13) ^ (^'""i^^r.vr. = h' atyeU with y„ > 



m=l 



where /i' is Holder continuous in the closure of U H {y : yn > 0} and 



(e+ < ^(a;)V-u, Vn >) 

Thus we see that (a^"*) as a matrix is Holder continuous and positive 
definite at each point. Consequently, the linear equations corresponding to 
(3.13) can be solved and vl'y^y^ is expressible in terms of functions which 
are Holder continuous in the closure of [7 R {y„ > 0}. Therefore u € C^'" 
upto yn = 0. Using the definition of u, we then obtain that u is C^'" in the 
closure of O fl (7(C/). Interior estimates are similar so we conclude that u is 
(j2,a £qj. gome a > in the closure of O n S(0, 2). 

We can now use Schauder type arguments as in [ADN] for equations and 

[ADNl] for systems (see also [Li] for quasilinear equations) to bootstrap 
and eventually deduce that u is infinitely differentiable in the closure of 
^ n 5(0.2). The proof of (2.7) is now complete. 



Remark: We emphasize that in order to satisfy the hypotheses in [Tl] (i.e 
to get an even Lipschitz extension of A) , it was important that ci;tn — 

for 
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k < n. This is precisely why we chose our coordinates using the distance 
function. 
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